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1. Introduction 

In the 1700's, Immanuel Kant and the Marquis de Laplace proposed that 
the solar system collapsed from a gaseous medium of roughly uniform den- 
sity ([79], [105]). A flattened gaseous disk - the protosolar nebula - formed 
out of this cloud. The Sun contracted out of material at the center of the 
disk, while the planets condensed in the outer portions. Despite its sim- 
plicity, this model suffered from the angular momentum problem inherent 
to star formation: a cloud of gas and dust with the diameter of the solar 
system and the mass of the Sun has too much angular momentum to col- 
lapse to the Sun's present size. It took 150 years to solve this problem. C. 
von Weiszacker worked out the basic physics of a viscous accretion disk, 
building on previous realizations that a turbulent viscosity could move ma- 
terial inwards and angular momentum outwards through the protosolar 
nebula ([178], [179]). With Liist's steady-state solution ([119]), Kant's neb- 
ular hypothesis no longer suffered from the angular momentum problem 
and became the leading model for solar system formation. 

Disk accretion is now known to power a wide range of astronomical 
sources. Starting with Kuiper's landmark study of [3 Lyr, viscous accretion 
disks have become central to our understanding of interacting binaries, 
where one star fills its inner Lagrangian surface and transfers matter into 
a disk surrounding a smaller companion ([103]). A disk surrounding a su- 
permassive black hole is widely believed to drive the phenomena observed 
in active galactic nuclei (e.g., [175]). Accretion disks remain an important 
issue in star and planet formation, as described throughout this volume. 

One striking feature of recent observations is the common phenomena 
observed in different accreting systems. All disks vary in brightness. These 
fluctuations range from rapid flickering on time scales as short as a few mil- 
liseconds up to long-lived eruptions with durations of decades or centuries. 
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Most accreting systems lose mass; practically all disks that lose mass eject 
material in a well-collimated jet ([115]). These common phenomena occur 
in systems with luminosities ranging from much less than 1 Lq in some 
interacting binaries up to 10^^ Lq in active galactic nuclei. 
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Figure 1 - Light curves for SS Cyg, a cataclysmic binary; CI Cyg, a symbiotic binary; 
and V1057 Cyg, an FU Ori variable star. Despite a different size scale in each system, the 
light curves have common features that indicate a similarity in their underlying physics. 



As an example of these common phenomena, Figure 1 shows light curves 
for SS Cyg - a cataclysmic binary with an orbital period of Port = 6.6 hr 
([30]); CI Cyg - a symbiotic binary with Po^ft = 855 days ([81], [16]); and 
V1057 Cyg - an FU Ori variable that is apparently a single star ([84]). 
The 3-5 mag eruptions of these systems occur when the accretion rate 
through the disk increases. Aside from the deep eclipses in CI Cyg, the 
three light curves have common features, with a rapid rise, a brief interval at 
maximum, and a long decay. All three systems resemble A- or F-type stars 
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at maximum light. This spectrum cools as the brightness fades. Substantial 
mass loss is associated with each system; the mass loss probably increases 
with the rise in brightness and decreases as the brightness declines. These 
similarities suggest that the same physical processes govern the evolution 
of disks ranging in size from 0.1 Rq to several tens of AU. 



TABLE 1. Types of disk systems 



System 


Distance (pc) 


Angular size (arcsec) 




BH + MS 


1000 


^io-« 


100 


NS + MS 


1000 


<io-^ 


100 


WD + MS 


100 


<io-3 


100+ 


MS + SG 


100 


^10-3 


10+ 


MS + RG 


1000 


^10-3 


10-100 


(3 Pic 


10+ 


1-10 


10-100 


Young star 


100+ 


1-2 


10-20 


AGN 


10* 


<>io-' 


1-10 



The common features of accreting systems provide good tools to test 
physical models of disks. Table 1 summarizes the main types of disk- 
accreting systems known today and provides simple comparisons of their 
angular sizes and the number of resolution elements across the face of the 
disk, Nres, that can be derived from modern observations. The first set of 
systems includes (i) compact binaries containing a disk surrounding a black 
hole (BH), a neutron star (NS), or a white dwarf star (WD), and (ii) wider 
binaries with a disk surrounding a main sequence star. A low mass main 
sequence star usually feeds the disk in compact binaries; a more evolved 
subgiant or red giant star feeds the disk in the wider binaries. In both types 
of binaries, the angular size of the disk is small compared to the resolution 
possible with current ground-based or space-based observatories. However, 
eclipses of the disk by the secondary star can be used to divide the disk 
into many resolution elements and to map out the physical structure of the 
disk with a few simphfying assumptions ([70], [71]). The dusty disks in the 
second group of systems are large enough to image directly at optical and 
radio wavelengths, although few systems have been mapped in much detail 
([164], [77]). The final group of active galactic nuclei have only recently 
been mapped with VLBI techniques; Nres will probably increase as these 
techniques become more sophisticated ([68]). 
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2. Steady Disks 

Accretion disks work to transport mass radially inwards and angular mo- 
mentum radially outwards. To understand how a disk manages this feat, 
consider a thin ring with two adjacent annuli at distances ri and r2 from 
a central star (Figure 2). Material in these annuli orbits the central star 
with velocities, vi and V2- The velocity difference between the two annuli, 
vi — V2 > 0, produces a frictional force that attempts to equalize the two 
orbital velocities. The energy lost to friction heats the annuli; some disk 
material then moves inwards to conserve total energy. This inward mass 
motion increases the angular momentum of the ring; some disk material 
moves outwards to conserve angular momentum. Energy and angular mo- 
mentum conservation thus lead to an expansion of the ring in response to 
frictional heating. The ring eventually expands into a disk, which generates 
heat fand radiation) at a level set bv the accretion rate. 



Figure 2 - Schematic view of two adjacent annuli in a disk surrounding a compact star. 
Annulus '1' lies inside annulus '2' and orbits the star (filled circle) at a higher velocity, 

Vl > V2- 

2.1. ACCRETION LUMINOSITIES AND TEMPERATURES 

Figure 3 shows the standard disk geometry. In a steady disk, material drifts 

radially inward at a constant rate, M. For an infinite disk, the total lumi- 
nosity generated by accretion is Lace = GM^^M / R^, which is 
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Figure 3 - Schematic view of a disk that extends to the stellar photosphere. Disk material 
drifts radially inward until it reaches the boundary layer, where the rotational velocity 
slows to match the stellar rotational velocity. 



in familiar units. As material drifts inwards, half of this accretion energy is 
radiated by the disk, 

_ 1 

I-' disk — • (2) 

The other half of the accretion energy becomes kinetic energy of orbital mo- 
tion around the central star. In most cases, disk material drifts inwards until 
it reaches the stellar photosphere (Figure 3). At this point, disk material 
with an angular velocity of = (GM^/i?^)^/^ must slow down to match 
the stellar angular velocity, 17^. This transition occurs in the "boundary 
layer," a narrow ring with a radial size, 

Ru4 ~ hl/R^ < i?^, (3) 

where hu is the local scale height (/iw/i?* < 1; [138], [140], [141], [132], [96], 
[143], [95]). The energy lost in the boundary layer is the difference between 
the rotational energy per unit mass in the disk and the rotational energy 
per unit mass in the central star: 

Lbl = \La.c ( . (4) 

^ V ^disk ) 

In some accreting systems, the magnetic field of the central star is strong 
enough to truncate the disk at a radius, i?m > -R*- A rough estimate of the 
truncation radius is the Alfven radius, 
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Figure 4 - Schematic view of a disk truncated by a strong dipolar magnetic field. The 
inner radius of the disk lies close to the corotation radius, where the disk and the stel- 
lar photosphere have the same angular velocity. Material falls from the inner disk onto 
the star along the magnetic field lines and produces two circular rings surrounding the 
magnetic axis. 



where //j, is the magnetic dipole moment of the central star. More rigorous 
estimates of Rm depend on the magnetic field geometry ([113], [46], [114], 
[29], [101], [191], [4], [107], [182]). The disk then has a total luminosity 



In this geometry, disk material falls onto the star along magnetic field 
lines (Figure 4). The infalling gas shocks and produces two rings at stellar 
latitudes, ±6, if the magnetic and rotational axes are parallel ([46], [101]). 
These rings are compressed into tilted ellipses if the magnetic axis is inclined 
with respect to the rotational axis, as expected in most systems ([88], [122]). 
In both cases, the luminosity of one ring is 






(6) 
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To derive the temperature structure of the disk and the boundary layer 
or magnetic accretion ring, we again use conservation of energy. The disk 
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temperature assumes an approximate balance between the energy lost from 
blackbody radiation and the torque due to frictional heating from mass 
motion through a distance AR in a gravitational potential: 

GMjl f 4 - IT^) « 2 ■ aTj (8) 



R + AR 
For AR <C R, this expression is 



A /? 

GM^M—y « AttRARxtT^ (9) 
R 



which results in 



The exact tempcratiu^c of a steady disk depends on how the torque 
transports energy through the disk. Standard models assume that frictional 
heating vanishes at the stellar photosphere, which leads to an extra term 
in the energy equation (8) above ([120], [161]). The temperature is then 
([13]): 

T^{R)=T,J-^\ l-^/^ (11) 



i4 



where 



r„cc = 13000 K 



10-5 M|yr-iy \ R^ ) 



(12) 



This expression for the temperature vanishes at the stellar surface, because 
the frictional heating vanishes. The temperature reaches a maximum at 
RjR^ = 49/36 ([13]): 

Different boundary conditions at the stellar surface yield slightly different, 
non-analytic, temperature laws ([134], [135]). These results are usually close 
to the standard temperature law, equation (11). 

The temperature of the boundary layer depends on the mass accretion 
rate. At high M, the region is optically thick ([120], [174]). Energy gener- 
ated in the thin "dynamical boundary layer" (equation 3) should diffuse 
over a broader ring with a size comparable to the thermal scale height, h\^\ 
([132]). The scale height is smaller than the stellar radius in most cases. 
This "thermal boundary layer" is then much larger than the dynamical 
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boundary layer and has a temperature 4-5 times larger than the maximum 
disk temperature ([91]). For low M, the boundary layer is optically thin and 
cannot radiate energy efficiently ([138], [140]). The thermal scale height is 
then comparable to the stellar radius, and the boundary layer resembles a 
hot stellar corona ([138], [140]). The details of boundary layer structure are 
sensitive to properties of the stellar photosphere and the viscosity mecha- 
nism, but these general considerations hold for many physical settings ([18], 
[17], [136], [72], [73], [47], [48], [137], [76], [97], [98]). 

The temperature of magnetic accretion rings also depends on geometry. 
If the magnetic axis is aligned with the rotational axis, the rings have equal 
luminosities and temperatures. The temperature is also constant with az- 
imuth along each ring ([122]). The two rings have unequal luminosities when 
the magnetic axis is tilted with respect to the rotational axis, as needed 
to produce the observed light variations in magnetic cataclysmic variables 
and T Tauri stars ([36], [88]). The luminosity also varies with azimuth along 
each ring ([122]). In the aligned case, the infalling matter shocks above the 
stellar surface and forms identical optically thick accretion columns that ra- 
diate in the Balmer and Paschen continua ([104], [26]). The surface area of 
this emission is comparable to the surface area of the boundary layer, so the 
temperature is still 3-5 times larger than the maximum disk temperature 
([26]). The accretion columns are probably quite different in non-aligned 
cases, but their structure has not been addressed. 



TABLE 2. Accretion and nuclear energies 



Object 




(erj 




enuc (erg g ^) 


T Tauri star 


5 


X 


10" 


5 X 10" 


main sequence star 


2 


X 




4 X 10^^ 


white dwarf star 


1 


X 


IQi^ 


4 X 10^* 


neutron star 


1 


X 


1q20 


6 X 10" 



To judge the effectiveness of accretion as an energy source. Table 2 
compares the gravitational potential energy, eg = GM^^/Ri,, of various 
1 Mq stars with the nuclear energy generation rate of the dominant fu- 
sion source, enuc- If a star can accrete and burn material from a disk, then 
the ratio of the accretion luminosity to the nuclear luminosity is simply 
the ratio of these two quantities, r = Lacc/Lnuc = (g/^nuc- This expression 
ignores the difference in accretion and nuclear time scales, but it provides 
a simple guide to the importance of accretion as a long-term energy source. 
Table 2 shows that accretion is rarely the major energy source for main 
sequence stars and white dwarfs; nuclear energy will always overwhelm ac- 
cretion by factors of 40—2000 on long time scales. Accretion is an excellent 
energy source for pre-main sequence stars and neutron stars. It is roughly 
10 times more effective than deuterium burning in pre-main sequence stars 
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and nearly 20 times more productive than complete hydrogen burning in 
neutron stars. 

2.2. TURBULENT VISCOSITY AND DISK TIMESCALES 

The source of the frictional heating in accretion disks remains controver- 
sial. Ordinary molecular viscosity is too small to generate mass motion 
on a reasonable time scale, ~ days for disks in interacting binary systems 
and ~ years to decades for disks in pre- main sequence stars. The large 
shear between adjacent disk annuli suggests that disks might be unstable 
to turbulent motions, which has led to many turbulent viscosity mecha- 
nisms. Convective eddies, gravitational instabilities, internal shocks, mag- 
netic stresses, sound waves, spiral density waves, and tidal forces have all 
been popular turbulence mechanisms in the past three decades (see, for ex- 



ample, [120], [161], [25], [112], [176], [7], [60], [172], [168], [24], [177], [157], 
[144], [5]). 



Recent work has shown that magnetic stresses in a differentially rotat- 
ing disk inevitably lead to turbulence ([9], [8], [110]). The growth time and 
effectiveness of these magnetohydrodynamic mechanisms make them the 
current leading candidate for viscosity in most applications. How this tur- 
bulence leads to significant mass motion in a real accretion disk remains an 
unsolved problem. 

Shakura & Sunyacv side-stepped the basic uncertainties of viscosity 
mechanisms when they developed the popular "ct-disk" model ([161]; see 
also [179]). In this approach, the frictional heating in the disk is due to a 
turbulent viscosity, u = acgh, where is the sound speed, h is the local 
scale height of the gas, and a is a dimensionless constant. This concept is 
similar to the mixing length theory of convection, with a serving the role 
of the mixing length. In most applications, a ^ 1-10; a needs to exceed 
^ 10~^ to allow material to move inwards on a reasonable time scale. 

This viscosity definition orders the important time scales for the disk 
([120]). The shortest disk time scale is the dynamical (orbital) time scale, 
which increases radially outward: 



The thermal time scale measures the rate that energy diffuses through the 



The thermal time scale is intermediate between the dynamical time scale 
and the viscous time scale, which measures the rate matter diffuses through 
the disk. 




(14) 



disk. 




(15) 
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Figure 5 - Schematic evolution of a point-like mass enhancement, a blob, in a viscous 
accretion disk. A central star lies at the center of the disk with an outer edge indicated 
by the dashed circle. The disk initially has a smooth radial decrease in the effective tem- 
perature, T, and the surface density, E (panel A). Differential rotation between adjacent 
annuli smooths the blob into a ring on the orbital time scale, To, which produces 5- 
function increases in T and S (panel B). The thermal energy in the blob moves inwards 
and outwards on the thermal time scale, tt, which produces a gaussian-likc perturbation 
in the temperature (panel C). The mass motion on this time scale is small; E changes 
little on the thermal time scale. The mass perturbation spreads out through the disk on 
the viscous time scale, tv, as indicated in panel D. 



^ ^ 40d^ / «_yV« (16) 



a V 1 / 

Figure 5 illustrates the evolution of the disk on these time scales. The 
top panel shows a blob of material superimposed on the smooth density and 
temperature structure of a steady disk. In this example, the blob has a mass, 
6M, larger than the mass of an annulus, and a thermal energy content, 6E, 
larger than the local thermal energy of an annulus. Disk rotation broadens 
the blob into a narrow ring in one or two rotation periods, which produces a 
narrow spike in the radial distributions of surface density and temperature. 
The thermal energy of this ring moves inwards and outwards on the thermal 
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time scale. The temperature distribution broadens considerably, but the 
surface density hardly changes. The mass in the ring finally moves inwards 
and outwards on the viscous time scale. This evolution produces an extra 
increase in the temperature due to viscous dissipation. 

Table 3 compares numerical values for the three disk time scales. These 
results assume = 1 Mq and a = 10~^. For most interacting binaries, the 
disk radius is ~ 0.1-10 R©; the time scales range from a fraction of a day for 
tb up to several tens of days for ry. The time scales increase dramatically 
for the larger disks in pre-main sequence stars and active galactic nuclei. 
The viscous time scale at the edge of a typical solar system is comparable 
to the disk lifetime. 



TABLE 3. Disk time scales 



Disk radius 


h/r 


T(K) 


Td 


Tt 


TV 


0.1 R© 


0.02 


2 X 10* 


0.004 d 


0.04 d 


20 d 


1.0 


0.02 


5 X 10^ 


0.120 d 


1.00 d 


400 d 


1 AU 


0.04 


250 


1 yr 


4 yr 


1 X 10" yr 


10 


0.06 


60 


32 yr 


105 yr 


1 X 10°^ yr 


100 


0.08 


25 


1000 yr 


2500 yr 


2 X 10® yr 



2.3. DISK ENERGY DISTRIBUTIONS 

Figure 6 shows the broadband spectral energy distribution of a steady ac- 
cretion disk with a boundary layer. The model assumes that the disk is 
optically thick and that each annulus radiates as a blackbody ([120], [13], 
[174], [158]). There are two main features in this spectrum: (i) radiation 
from the boundary layer at short wavelengths, and (ii) disk radiation at 
longer wavelengths. If the disk has a large ratio of outer radius to inner 
radius, Rout/Rin ^ !> the disk spectrum follows 

XFx oc < A < Xout (17) 

where XinTd^Rin) ~ XoutTd{Rout) ~ 0.36 ([120]). The disk spectrum follows 
the Wien tail of a hot blackbody, Td{Rin), at short wavelengths and the 
Rayleigh-Jeans tail of a cool blackbody, TdiRout), at long wavelengths. 

Although Figure 6 assumes a blackbody disk, the general features of the 
model change little if the disk radiates as some type of stellar atmosphere. 
The most important modifications for pre-main sequence stars involve ir- 
radiation from the central star and the infalling envelope, as described by 
Beckwith in this volume. 
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log Wavelength (|am) log Wavelength (|am) 

Figure 6 - Energy distributions for accretion disks. The left panel shows the energy 
distribution of a disk (dashed line), a boundary layer (dot-dashed line), and the total 
energy output (solid line) for M = 2 x lO"'^ Mq yr~^ onto a solar-type star. The right 
panel adds the energy distribution (dotted line) for a star with L* « L^. 



Lynden-Bell & Pringle first applied these concepts to disks in interact- 
ing binaries and pre-main sequence stars. They added radiation from the 
central star to the spectrum, as indicated in the right panel of Figure 6. 
This model suggests that identifying disk radiation can be difficult when 
the stellar luminosity exceeds the accretion luminosity. Table 4 compares 
temperatures and luminosities of pre-main sequence stars and steady disks 
for two accretions rates, M = 10~^ and 10~^ M0yr~^. At low masses, the 
contrast between the disk and star is considerable unless the accretion rate 
is ^ 10~^ MQyr~^. This contrast decreases for more massive stars. At 10 
Mq, the star is hotter and as luminous as the disk, because massive young 
stars burn hydrogen before reaching the main sequence. Table 3 shows that 
accretion cannot be a dominant energy source for main sequence stars; the 
difficulty of seeing an accretion disk against the background of a young O 
or B star therefore is not surprising. 



3. Unstable Accretion Disks 

Despite the elegant simplicity of the steady-state model, all disks vary their 
energy output. These fluctuations range from small, 10%-20%, amplitude 
flickering on the local dynamical time scale ([23]) to large-scale eruptions. 



13 



TABLE 4. Accretion and stellar luminosities in pre-main sequence stars 



M^Mq) 




i*(L©) 


10"® 


M©yr 1 

Lace (L©) 


10"^ 


Moyr 1 

Lace (L©) 


0.1 


3100 


0.1 


1850 


3 


5850 


275 


0.4 


3600 


0.5 


1825 


7 


5800 


675 


1.0 


4250 


7.0 


1800 


12 


5700 


1200 


3.0 


8650 


80 


1700 


25 


5350 


2350 


6.0 


17000 


1000 


2100 


50 


6650 


5000 


10.0 


21000 


5000 


1850 


60 


5850 


6000 


20.0 


29000 


45000 


1600 


75 


5000 


7450 



~ 3-5 mag or more, that can last for several times the local viscous time 
scale (Figure 1). In well-studied dwarf novae, the outbursts occur in cyclical 
patterns with quasi-periods of 10 or more days ([30], [130]). 

Osaki and Pringlc et al. first noted the possibility for thermal instabil- 
ities in accretion disks ([139], [131]). In the standard picture, the structure 
of a steady disk is derived from balancing heat generated by viscosity with 
radiative cooling. However, the viscous energy generation and radiative 
losses are set by local disk parameters; the input M is an external param- 
eter. Cooling can balance heating only if the local physical quantities can 
adjust to the input M. If the local physical variables cannot adjust, cooling 
cannot balance heating and a thermal instability is likely. 

Figure 7 displays a simple illustration of a thermal instability in a single 
disk annulus. The solid curves indicate loci of stable disks, where cooling 
precisely balances heating. To the right of these lines, heating exceeds cool- 
ing because viscous energy generation - set by the local surface density S 
through M = z^S - exceeds the radiative losses set by the effective temper- 
ature, Td. Cooling exceeds heating to the left of the stability lines. 

To illustrate the evolution of the instability, assume the annulus has (S, 
To) at A' and an input accretion rate. Mi . This input accretion rate exceeds 
the stable accretion rate at A', Mg = Mj^f . The annulus evolves up the solid 
line on the viscous time scale to find a solution where Ms = Mi. Before 
reaching this point, the annulus arrives at Mb, the last stable solution for 
this opacity source. At B, the annulus needs a larger surface density to 
accommodate the larger M, and a larger to radiate away the extra heat 
generated by the larger M. To do so, the annulus evolves to C; it makes 
this transition at constant S because the thermal time scale is shorter than 
the viscous time scale. Once at C, the annulus radiates more energy than 
generated by viscous dissipation [Ms = Mc > M\) and evolves down the 
solid line towards a solution where Ms = ^i- The annulus moves towards 
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lower S on the viscous time scale until it reaches point D, the last stable 
solution on this branch of the equilibrium curve. At D, the annulus drops 
to A to try to find the smaller Td needed to accommodate a smaller M . 
Once at A, the annulus retraces the steps around this 'limit-cycle' curve as 
long as M = Ml. 




Figure 7 - Outline of a disk instability in a single annulus. The disk evolves along the 
solid lines ABCD for an input M (horizontal arrow). Heavy solid lines correspond to 
stable equilibrium solutions; a dashed line indicates the locus of unstable equilibria. 

The evolution of an instability in a complete accretion disk follows the 
simple illustration. In a real disk, most annuli are close to B when a single 
annulus makes the transition from B to C. The increased temperature of 
a single annulus transfers heat to neighboring annuli on the thermal time 
scale; these annuli then jump to the 'high state' and propagate the eruption 
to their neighbors. Figure 8 shows this evolution for five snapshots in the 
evolution of a dwarf nova accretion disk. Time increases upwards in the 
Figure. The bottom panel shows the temperature distribution of the disk 
at the onset of the eruption. The inner edge of the disk jumps to the high 
state first; annuli at larger radii progressively follow until all of the disk 
resides in the high state in the top panel. The disk maintains this state for 
several days, and then retraces its path back to the low state. 
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9.0 9.5 10.0 
log Radius (cm) 

Figure 8 - Evolution of a disk instability in a dwarf nova disk. The bottom panel com- 
pares the radial temperature distribution of a time-dependent disk (solid line) with the 
temperature distribution of a steady disk (dashed line). The eruption begins in the next 
panel with a sharp increase in the temperature at the inner edge of the disk. This increase 
creates a 'wave' of increased temperature that propagates radially outwards in the disk 
as time moves forward (up in the figure) . In the top panel, the entire disk has reached the 
hot state; the actual temperature distribution (solid line) is then close to the steady-state 
temperature distribution (dashed line) 



One feature of the disk instability picture is that the disk is never in 
a steady-state. The dashed lines in Figure 8 plot steady-state temperature 
distributions normalized to the temperature at the inner edge of the disk. 
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The actual temperature gradient is much flatter than the steady-state gra- 
dient in the low state (bottom panel). The M through the disk increases 
radially outwards and is below the input M at each point in the disk. The 
disk approaches but never reaches the steady-state gradient in the high 
state (top panel). 

Actual disk temperature distributions have been derived from eclipse 
light curves similar to those shown in Figure 1 for CI Cyg. K. Horne &; col- 
laborators ([70], [71], [125], [10]) have used maximum entropy techniques to 
recover T(i[R) from multi-wavelength light curves and spectroscopic data 
assuming (i) the disk is azimuthally symmetric and (ii) the brightness tem- 
perature is close to the local blackbody temperature. Their results indicate 
that quiescent disks are rarely close to steady-state, with temperature dis- 
tributions usually much flatter than T(i{R) oc R^^^^. Systems in eruption ~ 
dwarf novae and symbiotic stars at maximum light - more closely resemble, 
but never achieve, the steady-state temperature distribution with q = 3/4. 
These results generally agree with the model predictions in Figure 8. 

4. Disk Eruptions in Pre-Main Sequence Stars 

Most pre-main sequence stars vary in brightness. Irregular brightness vari- 
ations of ;^ 1 mag are a defining feature of T Tauri stars and many Herbig 
AeBe stars ([78], [61]). Recent studies show that many of these variations 
are due to dark spots rotating with the stellar photosphere ([66], [21], [22], 
[67]) or bright spots produced at the base of a magnetic accretion column 
([21], [88]). Small eruptions of 1-3 mag lasting several years occur in the 
EXors, a poorly-studied class which includes EX Lup and DR Tau ([64]). 
More spectacular 3-6 mag eruptions occur in the FU Orionis variables, also 
known as FUors ([63], [59]). 

G. Herbig first associated the eruptions of FUors with pre-main sequence 
stars ([62], [63]; see also [3], [94]). FU Ori, which lies at the apex of a fan- 
shaped nebula within the dark cloud B35, brightened by a factor of over one 
hundred in 100-200 days (Figure 9; [180], [181], [62], [63], [74], [162], [75]). 
Thirty years later, Welin discovered the 5 mag brightening of V1057 Cyg 
within an eccentric ring of reflection nebulosity (Figure 9; [185], [63], [84], 
[99]). Herbig later noted the similarity between the optical spectra of these 
two stars with spectra of V1515 Cyg, a faint variable star embedded in arc- 
shaped nebulosity ([61]). He collected archival photographic photometry 
and identified a slow rise from rupg « 15.5 in the late 1940s to rupg « 13.5 
in the late 1970s ([63]). This brightness increase continued until 1980, when 
the star experienced a dramatic decline and slow recovery ([100], [90]). 

Herbig's demonstration that FU Ori - and by analogy other FUors - is 
a pre-main sequence star is straightforward. FUors are clearly associated 
with dark clouds, having radial velocities indistinguishable from the cloud 
velocity. The optical spectrum, including the high lithium abundance, has 
similarities with spectra of T Tauri stars. The event statistics are plan- 
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sible for pre-main sequence stars. Finally, the rise in brightness is a real 
luminosity increase that is not a nova outburst, the main alternative. 
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Figure 9 - Blue light curves for three FU Ori vaxiables. 



Many new FUors have been discovered since V1515 Cyg. The catalogued 
number now stands at 11 ([59], Table 5). Most have been observed to rise 
3-5 mag in optical or near-IR brightness in less than one year, as indicated 
by the 'Y' in the outburst (OB) column of Table 5. V1515 Cyg is the 
only known example to require a decade to rise to visual maximum, but 
the historical light curves for some systems are poorly documented. A few 
objects have been called FUors based on common spectroscopic properties, 
which are described more completely below. Most recent FUors arc more 
intimately associated with the densest dark clouds than the first members 
of the class, which suggests that many eruptions might have been missed. 
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There are also 3 candidate FUors listed as a second group of objects 
in Table 5. These objects have some properties in common with FUors, 
but more data are needed to see if they have the other characteristics as 
well ([169], [159]). A few other objects have one or two properties of known 
FUors, such as an outburst or unusTial hght curve (e.g., [41], [69], [2], ([163], 
[192]). These are more probably EXors or Ae/Be stars than FUors. 

4.1. BASIC PROPERTIES OF FU ORIONIS OBJECTS 

FUors share a distinctive set of morphological, photometric, and spectro- 
scopic characteristics (Table 5). Three are known binary stars, L1551 IRS5, 
Z CMa, and RNO IB/IC ([102], [86], [11], [171], [156]). Most have delicate 
fan-shaped or comma-shaped reflection nebulae (RN) on optical and near- 
IR images ([50], [148], [127], [118]). Optical jets and HH objects are common 
([148], [38], [59]). Most appear associated with large-scale molecular out- 
flows ([189], [43], [42], [121], [190], [44], [117], and references therein). These 
features - together with the broad, blue-shifted Na I and H I absorption 
features described below - demonstrate that FUors drive powerful winds 
that interact with the surrounding medium ([12], [37], [188]). 



TABLE 5. Selected properties of FU Orionis objects 



Object 


OB 


Opt ST 


IR ST 


RN 


HH/Jet 


Wind 


Outflow 


Radio 


L1551 IRS 5 


? 


K 


M 


Y 


Y 


? 


Y 


Y 


FU Ori 


Y 


F-G 


M 


Y 


N 


Y 


N 


N 


Z CMa 


? 


F-G 


M 


Y 


Y 


Y 


Y 


Y 


BBW 76 


@ 


G 


M 


Y 


? 


Y 


N 


? 


V346 Nor 


Y 


? 


C 


Y 


Y 


? 


Y 


? 


Par 21 


? 


F 


M 


Y 


Y 


Y 


? 


? 


V1515 Cyg 


Y 


G 


M 


Y 


N 


Y 


Y 


N 


V1057 Cyg 


Y 


G 


M 


Y 


N 


Y 


Y 


Y 


V1735 Cyg 


Y 


? 


M 


Y 


N 


? 


Y 


Y 


RNO IB/IC 


Y 


G 


M 


Y 


N 


? 


Y 


? 


IC 430 


? 


G 


? 


Y 


Y 


Y 


N 


Y 


PP 13S 


Y 


? 


M 


? 


? 


? 


Y 


? 


Re 50 


? 


G 


? 


Y 


Y 


Y 


Y 


Y 



FUors also have very unusual spectroscopic characteristics. The opti- 
cally visible sources have F-G, and sometimes K-type, giant or supergiant 
spectra ([63]; Tabic 5, see also [58], [149], [153])). The optical reflection 
nebulae of several embedded FUors also show G-type absorption features 
([167], [166]). All FUors but Z CMa and V346 Nor have very deep CO ab- 
sorption bands on near-IR spectra at 1.6 /um and 2.3 /um ([126], [40], [31], 
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[160], [82], [19]). These features resemble the CO absorption bands observed 
in red giants and are much stronger than those observed in any otlier pre- 
main sequence star ([56], [57], [51]). CO absorption in Z CMa is weakened 
by dust emission from an embedded companion; V346 Nor has weak CO 
emission instead of CO absorption ([102], [170]; [153]). Many FUors also 
display strong water absorption features, which strengthens the evidence 
for a ~ 2000 K photosphere ([27], [160], [28]). 

All FUors show large excesses of radiation over normal G supcrgiants 
at both ultraviolet and infrared wavelengths. The near-IR excess is clearly 
photospheric in origin, because the CO and H2O absorption features are 
so intense. The UV excesses in Z CMa and FU Ori appear associated with 
an A- or F-type photosphere that is hotter than the G-type photosphere 
observed at longer wavelengths ([89]). In addition to significant far-IR and 
submm emission ([183], [184], [84]), many FUors are strong radio continuum 
sources at cm wavelengths ([35], [154], [155], [116], [156]). This emission is 
not photospheric in some FUors and may be produced in the outflow or 
the jet. 

Herbig first noted broad absorption lines on optical spectra of several 
FUors ([63]). Hartmann & Kenyon later confirmed large rotational veloc- 
ities of V sin i ^ 15-60 km s-^ for V1057 Cyg, V1515 Cyg, and FU Ori 
([55]). This property is now characteristic of the class; all FUors have broad 
optical or infrared absorption lines or both ([58], [165], [166], [82]). The rota- 
tional velocity further depends on wavelength. The near-IR CO lines in FU 
Ori and V1057 Cyg have significantly smaller rotational velocities than the 
optical lines ([56], [57]). The optical rotational velocity smoothly increases 
with decreasing wavelength in Z CMa and V1057 Cyg ([186], [187]). In FU 
Ori itself, a powerful wind masks weak absorption lines and makes it difli- 
cult to detect any variation of rotational velocity with wavelength should 
one exist. 

Many FUors display doubled absorption lines on optical and near-IR 
spectra ([55], [56], [57], [58], [165], [166]). The two absorption components 
in V1057 Cyg are separated by 30-40 km s~^; these features have much 
larger separations in Z CMa and FU Ori. The long-term stability of the 
doubled absorption lines indicates that the lines are not produced by two 
stellar components in a binary system ([87]). 

Most FUors also show strong evidence for mass loss, in addition to 
the HH objects and CO outflows described above. Practically every FUor 
displays deep, blueshifted absorption components on Ha and Na I D ([12], 
[37], [39], [187], [54]); weak blueshifted CO absorption might be present in 
FU Ori ([57]). The optical line proflles can also change on month to year 
time scales ([12], [37], [187]). The near-IR line profiles may also change 
([19]). Both V1057 Cyg and V1515 Cyg have very pronounced line profile 
changes and dips in their optical light curves. These fluctuations suggest 
dust formation in a variable wind, but the data are not very extensive ([90]). 

Finally, FUor eruptions must be repetitive ([63], [55]). The event statis- 
tics of known FUors suggest that a young star must undergo 10-20 FUor 
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eruptions before reaching the main sequence. This estimate may be a lower 
limit, because some outbursts have certainly been missed. Reipurth's pro- 
posal that FUor eruptions produce HH objects suggests a recurrence time 
scale of ^ 1000 yr, based on the identification of multiple bowshocks (ejec- 
tion events) with dynamical separations of 500-2000 yr ([145], [146], [53], 
[152], [6]). The discovery of giant HH flows with dynamical time scales of 
10^-10^ yr ([151]) allows 10-100 eruptions if eruptions actually power the 
outflow (see below). 

4.2. FU ORIONIS OBJECTS AS ACCRETION DISKS 

The observations described above place severe constraints on possible FUor 
outburst mechanisms. With a luminosity of a few hundred (Table 6), 
a typical FUor emits 10^^ to 10^^ erg during the course of a 10-100 yr 
eruption. Outbursts also recur every lO'^-lO''^ yr. This process produces an 
object resembling a rapidly rotating F-G supergiant in the optical and a 
more slowly rotating M giant in the infrared. It must be unique to young 
stars, because we have not observed a FUor event in an older star system. 

Accretion is the most plausible energy source for FUor eruptions ([55], 
[109]). By analogy with dwarf novae, FUor eruptions are the high state of 
a thermal instability in a disk surrounding a pre-main sequence star. This 
instability occurs if the accretion rate from the molecular cloud core into 
the disk lies between the stable accretion rates in the low state and the 
high state (Figure 7). The system cycles between these two states as long 
as the surrounding cloud can replenish the disk between outbursts. 



TABLE 6. Disk Properties of FU Orionis Objects 



Object 


Lboi (L©) 


M*M(10-^ M|yr-i) 


Ji,(R©) 


L1551 IRS 5 


30 


0.1 


0.4 


FU Ori 


220 


1.7 


1.2 


Z CMa 


400 


4.0 


1.6 


BBW 76 


200 


1.4 


1.1 


V346 Nor 


160 


1.0 


1.0 


Par 21 


200 


1.4 


1.1 


V1515 Cyg 


135 


0.8 


0.9 


V1057 Cyg 


400 


4.0 


1.6 


V1735 Cyg 


200 


1.3 


1.1 


RNO IB/IC 


750 


3.5 


1.5 



Table 6 lists the disk properties needed to power observed FUor lumi- 
nosities from accretion. To make these estimates, I assume a steady disk 
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with Ldisk = 157 Lq{M^M/R^) and T^ax ~ 6500 K {M^M/R^y/^. Com- 
bining these and setting T^nax = 6500 K for an F-G supergiant star yields 

and 

These expressions produce the results in Table 6 for known FUors, assuming 
a typical inclination angle, cos i = 1/2. The median inner radius for the 
FUor sample is 1.1-1.2 R©; the median mass accretion rate is M^M ~ 
1-2 X 10~^ Mg yr~^. For a mass-radius relation for pre-main sequence stars 

with ages of 10^'' yr, ^ 0.2-0.3 M© and M « 0.5-1 x lO"^ Mgyr"!. 

This very simple analysis indicates that a young star must accrete mate- 
rial at ~ 10~^ Mq yr^^ during an eruption to power the observed luminos- 
ity. For an adopted outburst duration of 100 yr, this accretion rate results 
in a total accreted mass of ~ 0.01 Mq per eruption. The total mass accreted 
during a FUor eruption must be replenished during a low state lasting ~ 
10^-10^ yr, which results in an infah rate of ~ lO^^-lO""^ Mgyr'^ This 
range is close to the infall rates envisioned - and in some cases observed - 
for typical cloud cores ([1], [128], [124]). 



4.3. OBSERVATIONAL TESTS OF DISK MODELS 

Steady accretion disk models successfully explain many observations of 
FUors ([55], [87], [59]). Theoretical models of disk instabilities produce 
eruptions that generally resemble FUor events, although comparisons with 
observations are still in their early stages ([109], [32], [33], [80], [14], [15], 
[34], [129], [173]). The next few paragraphs summarize these results; Herbig 
&; Pctrov present a different interpretation ([65], [133]; see also [106]). 

Figure 10 shows light curves for V1057 Cyg, one of the best-studied 
FUors. These light curves closely resemble those of other accreting systems 
and provide strong support for an accretion model ([87], [33], [84], [15], 
[34], [111], [173]). The amplitude of the decline is largest in the UV and 
decreases monotonically from 0.3 fim to 5 fim. This behavior is character- 
istic of accreting systems, which evolve towards cooler temperatures as the 
luminosity declines ([91], [15], [111], [30]). The spectral type variation also 
agrees with disk model predictions. If we assign stellar effective tempera- 
tures to the optical spectrum at maximum (A5 star) and at the current 
epoch (G5 star), the UBV decline indicates a source with a roughly con- 
stant radius. This requirement is easy to achieve with an accretion disk if 
the radius of the central star remains constant. Most stars cannot evolve 
at constant radius if their effective temperatures change. 
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Figure 10 - Light curves for V1057 Cyg. The symbols alternate between filled circles (at 

B, R, and K) and plus signs (V and J). The system varied erratically at B = 14 17 prior 
to outburst and then rose 5-6 mag in less than one year. The visual observations indicate 
a similar rise time, although few prc-outburst observations are available. No pre-outburst 
near-IR data exist; K-band observations began shortly after the outburst was detected. 
The amplitude of the decline clearly depends on wavelength, with larger amplitudes at 
shorter wavelengths (see [84], [162], [99]). 



The light curves of other FUors also generally agree with disk model 
predictions. The rise times of 200-400 days in FU Ori and V1057 Cyg are 
comparable to the thermal time scale, as expected for an eruption that be- 
gins at the inner edge of the disk (Tabic 3; [15]). The several decade rise 
in V1515 Cyg is too long for such an "insidc-out" instability, but could be 
caused by an outburst that begins in the outer part of the disk . These 
"outside-in" eruptions can result from perturbations of the disk by a com- 
panion star or planet, in addition to the disk instability mechanism de- 
scribed above ([33], [20], [15], [34]). In all systems, the decay times are 
comparable to the viscous time scale (Table 3). 
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log Wavelength (|Lim) 

Figure 11 - Spectral energy distribution for FU Ori. The observations (filled circles) 
indicate a clear excess over a G-type star (dashed line). A disk model that produces a 
G-type optical spectrum accounts for most of the data (solid line). The modest far-IR 
excess over the disk spectrum at ~ 100 fiui is probably produced by the surrounding 
nebula responsible for the fan-shaped reflection nebula. 

Disk models also successfully explain FUor SEDs ([55], [87], [84]). Fig- 
ure 11 shows a dereddened SED for FU Ori. The deep UV and near-IR 
absorption features indicate that the UV and IR excesses over a G-type su- 
pergiant spectrum are photospheric ([27]). If disks radiate like stars at the 
effective temperatures appropriate for the observed absorption features, the 
SEDs for FU Ori and other FUors require that the surface area of the emit- 
ting region increases with increasing wavelength. In particular, the surface 
area of the near-IR source must be 10-20 times larger than the surface area 
of the optical source; the 300 K material responsible for the 10 fim excess 
must have roughly 100 times the emitting area of the optical continuum 
region. A disk in which the temperature decreases radially outward (§2) 
naturally explains this observation. Disk models account for the 0.3-10 fxm 
SEDs of V1057 Cyg, V1515 Cyg, and FU Ori quite well (Figure 11; see also 
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[55], [1], [87], [84], [15], [173]). Several other FUors with modest reddening 
have similar SEDs ([84]). 

Recent interferometric observations have resolved FU Ori at near-IR 
wavelengths. The angular size at 2.2 /xm agrees with predictions of the 
accretion disk model ([123]). 

The disk model fails to explain 10-100 fim data of many FUors, although 
the model fits both FU Ori and BBW 76 from 1-100 fim ([84]). However, 
the decline of the 10-20 fim light of V1057 Cyg follows the optical light 
curve very closely ([84]). This behavior suggests that the mid-IR radiation 
is optical light absorbed and reradiated by a surrounding envelope. An 
infall rate of 1-5 x 10~^ yr^^ produces an optically thick envelope that 
can reprocess optical light from the inner disk and account for the 10-100 
lim SED ([84], [15]). This rate is sufficient to replenish the disk in 1000 yr, 
which allows recurrent FUor eruptions in this system. The observational 
evidence for infall at comparable rates in two FUors, Z CMa ([108]) and 
L1551 IRS5 ([128], [124]), lends support to this interpretation. 

Finally, an accretion disk also naturally explains the optical and IR 
line profiles observed in some FUors ([56], [57], [87], [58], [186], [187]). The 
gradual decrease in the rotational velocity with increasing wavelength oc- 
curs because the longer wavelength emission is produced in more slowly- 
orbiting material at larger disk radii than the more rapidly moving inner 
disk material responsible for the short wavelength emission. At a given disk 
temperature, a larger fraction of the disk surface rotates at high line-of-sight 
velocities and the lines appear doubled ([87], [28], [137]). 

4.4. THE IMPORTANCE OF FU ORI ERUPTIONS 

Despite their relative rarity, FUors arc an important part of pre-main se- 
quence stellar evolution. FUor eruptions add a significant fraction of a 
stellar mass to the central star. FUor eruptions also eject a large amount of 
mass into the surrounding cloud. FUor eruptions may be frequent enough 
to provide nearly all of the mass of a typical low mass star. If this hypothe- 
sis is correct, FUor eruptions can power most molecular outflows and may 
be a dominant source of cloud heating in molecular cloud cores. This view 
is controversial, but can be tested with observations as outlined below. 

Although the statistics are still crude, young stars clearly accrete much 
of a stellar mass in FUor events. A typical eruption adds ~ 0.01 Mq to the 
central star; a rate of 10-20 eruptions per pre-main sequence star leads to 
a total accreted mass of ~ 0.1-0.2 Mq. These eruptions eject considerable 
amounts of mass and momentum into the surrounding cloud. For example, 
FU Ori itself has lost material with a momentum of ~ 0.3 Mq km s"-"^ since 
its eruption began. With 10-20 such eruptions, FU Ori could supply a small 
fraction of the momentum of a typical molecular outflow, {Mv)o ~ 1-100 
Mq km s-i ([45]). 

Two conditions must be satisfied for a young star to accrete all of its 
mass in FUor events: 1. The disk must undergo regular thermal instabilities 



25 



and cycle between the low and high states on a reasonable time scale. 2. 
FUors must be young, ^ a few x 10^ yr, to allow the surrounding cloud to 
fuel the recurring instabilities. 

As a group, FUors clearly are young objects. FUors have many more 
properties in common with young class I protostars than older optically 
visible T Tauri stars. FUors have distinctive reflection nebulae and large 
far-IR excesses; most arc associated with jets, HH objects, and molecular 
outflows. These properties are common in class I sources and rare in T 
Tauri stars ([145], [147], [49]). The typical age of a FUor is therefore ~ a 
few X 10^ yr ([85], [184], [84]). 

The conditions needed for the disk to cycle between the high and low 
states arc probably satisfied in many prc-main sequence stars. For a typical 
infall rate of 1-10 x 10~^ M0yr~^, theoretical models indicate that the 
disk spends most of its time in a low state where the mass accretion rate 
through the disk is much lower than the infall rate ([14], [15], [111])- The 
disk requires only ~ 10^-10^ yr to accumulate the ~ 0.01 Mq needed to 
power a FUor eruption. The disk will continue to cycle between the low 
and high states as long as the cloud supplies mass to the disk. 

If young stars accrete most of their mass in FUor events, it follows that 
FUors can power molecular outflows. In most outflow models, the mass 
ejected by the central star is ~ 10%-30% of the accreted mass. For a wind 
velocity of ~ 200 km s~^, the momentum in the wind is 

{Mv)u, ~ 20 (M^/Mq) M0 km s'^ (20) 

for a typical wind velocity of 200 km s~^. The observed range of outflow 
momenta - {Mv)o = 1-100 Mq km ([45]) - requires stellar masses of 
~ 0.1-5 Mq if the flows conserve momentum. The winds of FUors can 
power this outflow if a young star accretes most of its mass in FUor events. 

These simple estimates establish that FUors are important events if 
they recur on time scales of ~ 10^ yr during the protostellar phase of 
pre-main sequence stellar evolution. In this picture, the fraction of time in 
the FUor state is simply the ratio of the infall rate to the FUor accretion 
rate: / ~ Mi/Mpu- FUor disk models require Mpu 10~^ MQyr~^, so 
/ ^ 0.05 for a typical infall rate of Mj ~ 5 x 10~^ Mq yr~^. This prediction 
currently agrees with the statistics of FUors among protostars: 1 out of ~ 20 
protostars in the Taurus dark cloud contains a FUor (L1551 IRS5; [85]) and 
5 out of 17 protostars with HH objects contains a central star that resembles 
known FUors spectroscopically ([l50]; see also [52]). Sensitive photometric 
and spectroscopic surveys of nearby molecular clouds can improve these 
statistics considerably. If the FUor frequency among protostars turns out to 
be more than a few per cent, then FUors may represent the main accretion 
phase of early stellar evolution. 

I thank Janet Mattel for the SS Cyg data and John Cannizzo for the 
disk instability models. 
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